Quantum Phase Transition and Dynamically Enhanced Symmetry in 
Quadruple Quantum Dot System 
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We propose a system of four quantum dots designed to study the competition between three types 
of interactions; Heisenberg, Kondo and Ising. We find a rich phase diagram containing two sharp 
features: a quantum phase transition (QPT) between charge-ordered and charge-liquid phases, and 
a dramatic resonance in the charge liquid visible in the conductance. The QPT is of the Kosterlitz- 
Thouless type with a discontinuous jump in the conductance at the transition. We connect the 
resonance phenomenon with the degeneracy of three levels in the isolated quadruple dot and argue 
that this leads to a Kondo-like dynamical enhancement of symmetry from U{l)xZ2 to C/(l)xC/(l). 

PACS numbers: 73.21. La, 05.30.Rt, 72.10.Fk, 73.23.Hk 



Strong electronic correlations create a variety of inter- 
esting phenomena including quantum phase transitions 
[H, emergence of new symmetries [2|,non- Fermi-liquid 
behavior and Kondo resonances Q. It is likely that 

new, yet undiscovered, phenomena can arise from unex- 
plored competing interactions. Today, quantum dots pro- 
vide controlled and tunable experimental quantum sys- 
tems to study strong correlation effects. Further, un- 
like most materials, quantum dots can be modeled using 
impurity models that can be treated theoretically much 
more easily. Single quantum dots have been studied ex- 
tensively, both theoretically and experimentally, which 
has led to a firm understanding of their Kondo physics 
[1, ■ More recently, the focus has shifted to multiple 
quantum dot systems where a richer variety of quantum 
phenomena become accessible Q • These include non- 
Fermi liquids dynamical enhancement of symmetry 
0, and quantum phase transitions lol- 1^. 

In this work we propose a quadruple quantum dot sys- 
tem, that is experimentally realizable, in which three 
competing interactions determine the low temperature 
physics: (1) Kondo-like coupling of each dot with its lead, 
(2) Heisenberg coupling between the dots, and (3) Ising 
coupling between the dots. Thus, there are two dimen- 
sionless parameters with which to tune the competition. 
The pairwise competing interactions, Kondo/Heisenberg 
and Kondo/Ising, have both been studied previously. 
The two impurity Kondo model with a Heisenberg inter- 
action between the impurities shows an impurity QPT 
from separate Kondo screening of the two spins at 
small exchange to a local spin singlet (LSS) phase at 
larg e exchange. This has received extensive theoretical 
[id I3I 1^ and experimental ?| attention. The compe- 
tition between Kondo and Ising couplings has also been 
studied theoretically for two impurities [ll|, including 
in the quantum dot context ll|, ll2|; however, no experi- 
mentally possible realization of this competition has been 
proposed to date. 

Our system consists of four quantum dots and four 
leads, as shown in Fig.[lja), with two polarized (spinless) 



electrons on the four dots. We find that the system has a 
rich phase diagram, Fig.[TJb), in terms of the strength of 
the Heisenberg interaction controlled by t and the Ising 
interaction controlled by U' . In the absence of the Ising 
interaction we start in the local spin singlet (LSS) phase. 
Upon increasing the Ising strength, we find that the sys- 
tem first evolves continuously to a new Kondo-type state 
with a novel C^(l) x U{1) strong coupling fixed point. 
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FIG. 1: (color online) (a) Quadruple-dot system. U and U' 
are electrostatic interactions while t and F involve electron 
tunnehng. (b) Ground state phase diagram as a function 
of the Ising/Kondo tuning f/'/F and the Heisenberg/Kondo 
tuning t/T. The two distinct phases — charge ordered (CO) 
and charge liquid — are separated by a KT quantum phase 
transition [blue crosses (numerical) and green dashed hue 
(schematic)]. Several cross-overs lie within the charge- liquid 
phase. Red stars mark the level crossing where the U{l)x (7(1) 
state is found (numerical). The charge Kondo region lies be- 
tween the red and blue lines. "LSS" denotes the local spin sin- 
glet state (Heisenberg coupling dominates), while when both 
Heisenberg and Ising couplings are weak, the system consists 
of inidividually screened Kondo states on the left and right. 
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Then there is a crossover to a SU{2) charge Kendo state. 
Finally, an additional small increase in U' causes a QPT 
of the Kosterlitz-Thouless (KT) type to a charge ordered 
state (CO) (as in Refs. [lll[l3 ) consisting of an unscreened 
doubly degenerate ground state [isj . 

Model — The quantum dots in Fig.[IJa) are capacitively 
coupled in two ways: U is the vertical interaction (be- 
tween L+ and L— ; i?+ and R—) and U' is along the 
diagonal (between i+ and i?— ; L— and R+). Along the 
horizontal, there is no capacitive coupling but there is 
direct tunneling t (between L-\- and -R-I-; L— and R—). 
Each dot couples to a conduction lead through F — nV^p 
where p is the density of states of the leads at the Fermi 
energy. The whole system is spinless. We consider only 
the regime in which the four dots contain 2 electrons. 

The system Hamiltonian is H = i?iead + -f^imp + ^^coup , 
where -ffiead = J2i.s,k <^kclg,,Ctsk describes the four con- 
duction leads (i = L,R; s = -I-,—), and i?coup = 
s fc('^Lfc'^is+H.c.) describes the coupling of the leads 
to the dots which produces the Kondo interaction. iJimp 
is the Anderson-type Hamiltonian 



imp 



i=L,Rs=+,- i=L,R 
+ U'iflL+flR- + nL-nR+) + * X! ^^Ls^-Rs 



(1) 



it c 



We take J7 3> C/' so that there is one electron on the left, 
and one on the right. 

We can reformulate ffimp as an exchange Hamiltonian 
by noticing that the right-hand (left-hand) sites form a 
pseudo-spin: Si = J^s s' dl^(Tss'dsi/2. When t <Si U , the 
effective Hamiltonian for the quantum dots is 



off 



Jh5'l • Sr — JzS^Sr 



(2) 



where Jh At^/{U - U' /2) and ~ 2U' . Thus t con- 
trols the strength of the Heisenberg interaction among 
the dots, and U' controls the Ising coupling. The eigen- 
states of the impurity site are the usual (pseudo)spin sin- 
glet and triplet states, \S), \++), |~-), and |rO). 

Two limits of our model have been studied previously. 
First, for U' = 0, it becomes the well-known two im- 
purity Kondo model [13, M\- If direct charge transfer 
is totally suppressed, a QPT occurs between a Kondo 
screened state (in which the impurities fluctuate between 
all four states, singlet and triplet) and a local spin singlet 
(LSS) 13, When direct tunneling is introduced, the 
QPT is replaced by a smooth crossover . Second, when 
t = 0, the model has [Til, [13] a Kosterlitz-Thouless (KT) 
type QPT between the Kondo screened phase at small 
U' and a charge ordered (CO) phase at large U' . The 
CO phase has an unscreened doubly degenerate ground 
state corresponding to and | ). 

We solve the model ^ exactly by using finite-temp- 
erature world line quantum Monte Carlo (QMC) simula- 
tion with directed loop updates [13, [13 • We study the 
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FIG. 2: (color online) (a) Local charge susceptibility as a func- 
tion of temperature. The power-law behavior of the top three 
curves indicates the CO phase. The peak and low-T constant 
in the lowest curves indicate the LSS state. The low-T sat- 
uration of the middle curves is due to Kondo-like screening, 
(b) Pesudospin-pseudospin correlation as a function of C/' for 
different /3. Inset: Zoom near the crossing point. The cross- 
ing of the singlet and doublet levels occurs at U[^q = 0.142. 
corresponding to level crossing. ([/ = 3, F = 0.2, and t = 0.3.) 



LSS state in the absence of 



M'^jU > T^'^ where T^'^ is the Kondo 



regime in which there is 
Ising coupling, / l/ ^ vvncic 
temperature of the left or right pseudospin individually. 
Taking the leads to have a symmetric constant density 
of states, p = 1/2D, with half-band-width D = 2, we 
focus on the case [/ = 3, F = 0.2, and t = 0.3. ,3 = l/T 
is the inverse temperature. As U' is varied [a horizontal 
scan in Fig.[T](b)], the gate potential is chosen such that 

= ^{U + U')/2, placing the dots right at the midpoint 
of the two electron regime. 

Thermodynamics — As a first step toward distinguish- 
ing the different phases, we look at the local charge sus- 
ceptibility xi°'' = Io{MT)MO)}dT, where A = nL+ + 
nfl-i- — riL- — n/f_ and n.i s is the charge density of the 
dot labeled i, s. This could be measured experimentally 
using recently developed single charge sensing techniques 
[l3] . Fig.[2Ua) shows Xc°'^ as a function of temperature for 
different values of U'. The curves show three types of be- 
havior. First, for small Ising coupling ([/' < 0.11), Xc°'^ is 
roughly constant at low T and has a peak at higher tem- 
perature. This is the LSS phase. The value of T at which 
xj.°'^ peaks decreases as the energy spacing between the 
singlet \S) and doublet, {|-l~l-), | — )}, decreases. Second, 
at the other extreme, for large Ising coupling {If > 0.15), 
xj,°'^ behaves as 1/T down to our lowest T. This is a clear 
signature of the CO phase in which the two charge states 

|-|— 1-) and I ) are degenerate. The tunneling t does not 

lead to any relevant operator which might destroy the CO 
phase [lij : thus, although the phase boundary is a func- 
tion of t, the essential nature of the KT QPT [Fig.[T];b)] 
is not affected. Third, for intermediate values of U' , Xc°'^ 
becomes large and then either decreases slightly at our 
lowest T or saturates. This behavior can be produced 
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by either a near degeneracy between the singlet and dou- 
blet states or by charge Kondo screening of the doublet 

{|++), I )}. As we will see from the conductance data 

below, the QPT to the CO phase occurs at a value J/^rp 
between 0.146 and 0.15. 

To extract the position of the level crossing between 
\S) and {! + +), I )}, we calculate the pseudospin cor- 
relation function (SlSfi) as a function of U' for different 
T [Fig.Hb)], where 5f = {n,+ - nj_)/2. For U' = 0, the 
ground state is the LSS so that {Sf^S^Pj) ~ —0.2 is close to 
— 1/4. On the other hand, for large C/', in the CO phase, 
{^l^r) is positive and approaches 1/4. (The charge fluc- 
tuations due to tunneling to the leads causes the values 
to differ slightly from ±1/4.) The crossing point of the 
curves for different temperatures gives the position of the 
(renormalized) level crossing. The inset shows that it 
occurs at {SIS^) « 1/12, which is consistent with the 
isolated-dots limit. The position of the level crossing is, 
then, [/lq « 0.142; note that this does not coincide with 
the QPT to the CO phase (0.146 < t/^T < 0.15). 

Conductance — Conductance is a crucial observable ex- 
perimentally. However, QMC is only able to provide nu- 
merical data for the imaginary time Green function at 
discrete Matsubara frequencies — the conductance cannot 
be directly calculated. The zero bias conductance for an 
impurity model can be obtained [l9| by extrapolating 
to zero frequency: G ~ lim^j^^^o g(ia;„) where g{iuJn) = 

dTCOs{ujnT){Pa;{T)Py{0))ujn/h and Py = J2y'>yf^V'- 
We have recently shown that this method works very well 
for Anderson-type iinpurity models in the Kondo region 
at low temperature [20|. 

We use this method [ll] to find the conductance be- 
tween the left and right leads as a function of U' for 
different T; the results are shown in Fig.|31 For U' 
small ([/' < 0.1), the conductance is small because the 
phase shift is nearly zero in the LSS state [l3|- For U' 
large {W > 0.15), the conductance is also small and ap- 
proaches zero as [/' — > oo, consistent with the argument 
in Ref. 11. At intermediate values of U', there is a strik- 
ingly sharp conductance peak near the value of U' where 
the level crossing occurs. Here, the conductance increases 
as T decreases and approaches the unitary limit 2e^//i as 
T — )• 0. The position of the conductance peak approaches 
the level crossing U' = 0.142 at low temperature [15|- 
association with the level crossing suggests that this peak 
comes from fluctuations produced by the degeneracy of 
|5) and{| + +),|— )}. 

A sharp jump appears after the peak: notice that the 
conductance at U' — 0.146 increases at lower tempera- 
ture while that at 0.15 decreases [see panels (b) and (c) 
for clarity]. The latter behavior is the signature of the 
CO phase, while the former suggests a Kondo-like phase, 
namely the dynamic screening of the | )} dou- 

blet. Thus, this sharp jump is associated with the KT 
QPT from the screened to the CO phase [llj, which oc- 
curs between U' = 0.146 and 0.15. 
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FIG. 3: (color online) (a) Zero bias conductance as a function 
of U' for two values of /3. Inset: Zoom on the peak caused by 
the U{l)x U{1) ground state. The T = expectation from the 
effective theory near the level crossing is indicated schemat- 
ically by the black solid line; the two points of discontinuity 
(the level crossing and the KT QPT) are marked by dashed 
lines. (b),(c) Conductance as a function of temperature for 
U' = 0.146 and 0.15, respectively; the opposite trend in these 
two curves shows that they are on opposite sides of the QPT. 

Effective theory near the level crossing — To gain in- 
sight into the conductance peak, we develop an effective 
theory near the level crossing. Using T/U as a small 
parameter, we make a Schrieffer- Wolff transformation to 
integrate out |rO); to include tunneling, processes of or- 
der Tt/U^ must be included 2l|. Higher-order terms in 
T/U arc neglected. In the leads, only the combinations 



Cq need be considered as these are the lo- 



cations to which the dots couple. The resulting effective 
Kondo Hamiltonian reads 

^Kondo = A{MlSl+MlS\) + 2JlMlSl 

+ jf{M]^S'^} + M^}S^l)+2J^^M]^S^^. (3) 

The operators M act on the dots, 

= V2{\++){S\ T |^)(- 
Ml = I ++X++I - I— >(- 
M'^ = | ++>(++| + I— X- 

while the operators S act on the lead sites. 



-|) = (Afi/")t 

-I 

-|-2|^)(^|, (4) 



-0,1- 



.Co,i-) 



i=L,R 



(5) 
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For t/U <^ 1 and particle-hole symmetry, ~ ~ 
W^/{U + U') and Jii ~ Jfi ~ 8F2t/(f7 + U'f. 

Renormalization effects in -ffKondo '^^^ analyzed us- 
ing poor man's scaling 22], yielding the scaling equations 



dj\/d\nD = 

dfl/dhiD = 

dJl/dlnD ^ 

dJl^/d\nD = 



-2p{JiJl + 3j'2jI') 

-2p[(Ji)2 + (Jii)2] 



(6) 



Numerical solution of these equations reveals that at a 
certain value of D, all the coupling constants simulta- 
neously diverge. This defines the problem's characteris- 
tic energy scale Dq, which can be considered the Kondo 
temperature. The coupling constants have a fixed ratio 
as they diverge: limo-^Do J]_ '■ Ji: ■ Jl '■ Jl^ ^ \/2 : \/2 : 1 : 
1, suggesting dynamical symmetry enhancement of the 
ground state. 

Symmetry analysis — The six S operators form an 
SO{A) algebra 1^. However, the six M operators do not; 
rather they form part of an SU{3) algebra — the missing 

operators are | + +)( | and | )(++! U^]- Since 

^Kondo the product of two objects which generate dif- 
ferent algebras, the symmetry of the system must be a 
subgroup of both 50(4) and SU{3). To study the com- 
plete symmetry group of both the bare and fixed-point 
Hamiltonians, consider the total z-component of pseu- 
dospins type I and II 



5'z,Tot = Ml 



where sl^j^ is defined by replacing cq with Ck in Eq. ([5|) . 

One can check that [•S'^^t^j, i?iead + -ffKondol = 0' which 
gives a (pseudo)spin C/(l) symmetry for the bare Hamil- 
tonian. The bare Hamiltonian also commutes with inter- 
changing L and R or interchanging + and — . Thus, the 
symmetry of the bare Hamiltonian is U{l)sx .^2,l_r x 
[an irrelevant charge U{1) is ignored]. 

At the fixed point, lim^j-j/j,-, J]_/ J± — >■ 1 implies that 
both Sl rp^^. and S'^^rp^j commute with the Hamiltonian. 
Thus there is an additional U{1) symmetry. Note that 
exp{i9S]^'Y^^) generates the L ^ R transformation for 
9 = TT. Therefore, the ^2,l_r symmetry of the bare 
Hamiltonian is dynamically enhanced to a U{1) symme- 
try. The complete symmetry group at the fixed point 

(ground state) is U{l)s x U{1)lr x ^2,^ , where the 

2^2, H symmetry is irrelevant for the Kondo physics. 

Experimental accessibility — Though our analysis above 
is for a highly symmetric model, possible experimental 
observation is greatly aided by the fact that this symme- 
try is not essential. For our scenario, it is crucial that 

]-!--!-) and ] ) be degenerate, which can be achieved 

by fine tuning the gate voltages. When the levels are 
only approximately degenerate eL+ — £l- — e_R+ — e_R- 
(define Se as the detuning from the symmetric point). 



5e induces a relevant perturbation in the CO phase tl2|. 
However, a continuous but sharp crossover does still oc- 
cur in conductance as long as 5e <^ T^*^ and T > Se 
[12], where T^'-^ is the Kondo temperature at the level 
crossing. Other deviations from symmetry — those in- 
volving the tunneling between the dots or the coupling to 
the leads, for instance — merely change the coupling con- 
stants in the effective Hamiltonian. In both cases, the 
U{l)x U{1) strong-coupling fixed point remains stable. 
Finally, note that the transition also appears by tuning 
the experimentally more accessible t at fixed U' 15| [i.e. 
a vertical cut in Fig.[T][b)]. 

This work was supported in part by the U.S. NSF 
Grant No. DMR-0506953. 
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Supplementary Material for "Quantum Phase Transition and 
Dynamically Enhanced Symmetry in Quadruple Quantum Dot System" 

In this supplementary material we address (1) extracting the conductance, (2) the shape of the conductance peak, 
(3) more detail about the effective theory near the level crossing, and (4) an effective theory near the QPT. 



I. CONDUCTANCE FROM QUANTUM MONTE CARLO 



In this section, we show how to obtain the conductance from quantum Monte Carlo data by the extrapolation 
method; a longer description, as well as checks, is in Ref. The quadruple quantum dot system can be mapped 
to a one-dimensional infinite tight-binding chain as shown in Fig.|4l where a pseudospin is considered on each site 
(corresponding to + or — in Fig. 1 of the main text). The linear conductance is obtained [H, by extrapolating the 
conductance at the (imaginary) Matsubara frequencies, g{iujn) with a;„ = 27mr, to zero frequency, G = lim^^_>o 5(*w„), 



g(iw„) = / dTCOs(ijJnT){Px{T)Py{0)) 



(8) 



where Py is the sum of the electron charge density operators to the right of y, Py = J2y'>y ^y' ■ combinations 
of X and y can be used in Eq. ([8]) because the system is not a physical chain, but only effectively mapped to a chain. 
Notice that the current through the five bonds closest to the quantum dot sites (labeled LO and i?0) correspond to 
the physical current. Therefore, x and y must be chosen from among {LI, LO, i?0, i?2}. In addition, left-right 
symmetry reduces the number of independent combinations. In our calculation, we choose six cases for x and 
(LO,LO), (i?0,i?0), (L1,L0), {Ll,m), (Ll,i?l), and (LO,i?0). We carry out this extrapolation as in Refs.i and 
For the CO and LSS regions, the first several QMC data points [g{iijjn), n — 1, . . . ,m, where m = 4 to 6 depending 
on the case] can be fit to a quadratic polynomial, as in Fig. [5ja) (b) . For the Kondo regime, the first 14 QMC data 
points [g{iijjn), n = 1 ... 14] are fit to a series of rational polynomial functions [see Fig.[5l[c)(d) for examples]. The 
conductance is the average of the different extrapolation results, and the error bar corresponds to the maximum 



L3 L2 LI LO RO Rl R2 R3 
FIG. 4: (color online) The ID infinite tight-binding chain, where LO and _R0 mark the impurity (quantum dot) sites. 
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FIG. 5: (color online) Conductance at Matsubara frequencies (symbols) for quadruple dot system and the corresponding fits 
used to extrapolate to zero frequency (lines). In the calculation, the tunneling t — 0.3 and the coupling F — 0.2. The values of 
U' and /3 are (a) 0.2, 1000, (b) 0.11, 1000; (c) 0.137, 1000; (d) 0.141, 2000. 
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FIG. 6: (color online) Zero bais conductance as a function of U' for different /3 in the regime of the peak caused by the 
U{l)x U{1) ground state. The two vertical dashed lines show the level-crossing and QPT points. The black solid line gives the 
conductance at zero temperature schematically. 

spread. For all of the cases studied in this work, the extrapolation appears to be straight forward and reliable (for 
cases in other systems for which the extrapolation does not work well, see Ref.[ll). 



II. PEAK OF CONDUCTANCE NEAR RESONANCE 



Fig. [6] shows that the conductance peak approaches the level-crossing point as the temperature approaches zero, 
as expected from the T = theory, but that the peak is not at the level-crossing point at non-zero T. This shift 
in the peak position can be understood as follows. When the parameter U' increases toward the level-crossing point 
L^Lc ~ 0.f42, the singlet 15) is the ground state of the dots, and the level splitting A between \S) and the doublet 

{1++); I )} states decreases to zero at J/lc- Following the argument in Ref.[^ concerning the effect of a magnetic 

field on the usual quantum dot Kondo effect, we expect that the effect of the level splitting A at zero temperature 
should be similar to that of the temperature at zero level splitting, though G{A/Tk,T — 0) and G{T/Tk, A — 0) 
may have different universal forms. For finite A <C T, G can be approximated by G(T/Tk, A = 0) which saturates 
at 2e^//i at low temperature. On the other hand, at T = 0, G{T = 0, A) < 2e'^/h (because the LSS has a substantial 
weight in the ground state). Therefore, the conductance G{T/Tk, A) may show a peak at non-zero T, and does so for, 
e.g., {/' = O.f 35 in Fig. (5] Similarly, for fixed T ^ 0, G{A/Tk) may peak at non-zero A. In summary, although the 
conductance peaks at A for T = 0, at non-zero temperature, G{U') may show a peak at U' < 0.142 corresponding 
to finite level splitting A. 



III. EFFECTIVE THEORY NEAR THE LEVEL CROSSING 



The S operators defined in Eq. (5) of the main text satisfy the following commutation relations: 

[S\,Sl]=2Sl [SlSl] = ±Sl 

[S'lS'}] = 2Sl [S'^^,Sl] = ±S\ 

\S\ , S^l] = 23"^ [Si , S^] = ±Sl (9) 

[S'lSl] = 2Sf [Sf,Sl]^±Sl 



These relations generate the SO{A) algebra, so the six S operators form an SO{A) algebra. However, the six M 
operators do not form an SO{A) algebra. If the standard basis for the fundamental representation of SU{'i) is F^, 
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i = l,2,...,8[l], the six M operators can be written as the hnear combinations 
Ml = F3, Mf = (2/V3)f^8, 

Ml^{F^ + iFr,~F^ + iF^)/V2, Ml^{Fi~iFz~F^~iF^)/V2, (10) 
M:|i = (F4 + iFz + Fe - iFy)/^!, Af" = (F4 - ii^g + i^e + iF7)/V2 . 

The two missing operators | + +)( 1 and | )(+ + | can be written as 

! + +)( |-Fi+iF2, I ){++\^Fi-iF2. (11) 

Therefore, the six M operators combined with the two missing operators form an SU{?>) algebra. 

Up to order T /U , there is no direct process which leads to |++) \ ). The only path leading to | ++) <S=> | ) 

is via the singlet state as an intermediary: | + +) 4^ \S) ^ \ ). When higher-order terms (T/U)'^ are considered, 

the fom' electron hopping terms do produce a direct | + +) <S=> | ) process: 

^mp = ^ E (l + +>(— l4L-C,-L+4i^-Vi^++h-C-) ■ (^^^ 

k,k ,q,q 

Since is comprised of four electron operators, its naive scaling dimension is negative, suggesting that is 

irrelevant. To check this, consider the one-loop RG equations. The scaling equations for the system combining -ff^? 
with -ffxondo [Eq. (3) of the main text] consist of the original four equations in Eq. (6) of the main text plus one 
additional equation: 

= -2pJiA. (13) 



dlnD 



Solving these five equations numerically with Jl_ > 0, we find that the coupling A flows to for any initial value. 
Therefore, -ff^? is an irrelevant operator in the strong coupling phase, which confirms the naive scaling-dimension 
analysis. 

IV. EFFECTIVE THEORY NEAR THE QPT 

To study the physics near the KT quantum phase transition, we develop a low energy effective theory closely 
following Refs.0 and0. First, consider the effective Hamiltonian in the large U' limit without tunneling, and note 

that the energy of jiS*) is much higher than that of {| + +),| )}• Using T/U as a small parameter, we make a 

Schrieffer- Wolff transformation to integrate out \T0) and \S); higher-order terms are neglected. The resulting fixed 
point Hamiltonian is 



^CO - E ^ksact^Cksa + K ^ [flLa + flR, - ljc+ ^C^',^ . (14) 

k,s=L/R,a=+/~ kk'scr 

The ground state corresponds to a charge ordered state with two-fold degeneracy, {| + +), | )}: the two electrons 

are frozen in either the upper two quantum dots or the lower dots. The charge order can be screened by ifg? given 
in Eq. The model here is very similar to that for the orthogonality catastrophe in the x-ray edge problem 0,0], 
as pointed out in Refs.0 and@; we briefly summarize their argument here. 

The charge ordered state can be flipped by the operator / = |++) ( 1 c^i_Cj.'^^c^^_c^'^_,_. According to Hopfield's 

rule of thumb loj , the correlation function in the charge-ordered phase is given by 

/+(i)/(0))^^„ ^i"" where ^ (An,)' (15) 

i=L+,L-M+,R- 

is related to the change in occupation An,; of each dot. An^ can be expressed in terms of the conduction band phase 
shift d through the Friedel sum rule; therefore, the anomalous exponet can be related to 6, 

Y: (An,)^=4f^-l)' . (16) 



i=L + ,L-,R+M- 
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In a ID problem such as ours, the power-law decay of correlations in Eq. ([T5|) leads to the criteria 

irrelevant Charge Ordered state 
relevant Charge Kondo state 
marginal Critical Point 

because of the possible infrared (long time) singularity. Thus the criterion for the critical point separating the 
charge-ordered and charge-Kondo states is 

The phase shift S depends on the coefficient K of the potential scattering term in Hqq, while K itself depends on U, 
U', and r. Therefore, the phase shift is a function oi U, U' , and F, and 

<5,(c/,c/^T,r) = f (i-^) => c/^cT = /(f/,r). (18) 

So, when U > I7^ti '^^ have a > 2, H^fp is irrelevant, and the system is in the charge-ordered phase with a two-fold 
degenerate ground state. However, when U < ?7ktj < 2, H^f^ is relevant, and so this four electron hopping operator 

prodcues |-|--l-) \ ) which screens the charge-ordered state to form the charge-Kondo state. There is a KT-type 

quantum phase transition 0,0 at U' = U^^. 

In contrast to Refs.|^ and U we must consider the influence of the tunneling t on the QPT (i.e. direct tunneling 
between L+ and R+, and between L— and R—). To include tunneling, processes of order Tt/U^ need to be considered 
[lo| . The low energy effective Hamiltonian in the charge-ordered phase becomes 

k.s=L/R,a=+/- kk' sa 



Kt J2 ("^-^ + "^-^ " l)(4L^Cfc'fl^ -I- c+R^Cfc'i^) . (19) 

k,k\ a 

The contribution of the tunneling merely adds the possibility of potential scattering between the electrons in the left 

and right leads; note that the dependence on the filling of the four dots, either | + +) or | ), is the same. Thus, 

tunneling does not introduce a charge order flip process, but does contribute to the phase shift experienced by the 
lead electrons when the dot flips. We expect, then, that t does not affect the essential physics of the charge-ordered 
phase and QPT. We should check, however, whether the tunneling leads to a relevant process in addition to H^f^ 
which can screen the charge order. The lowest order screening terms induced by tunneling are 

A - I + +>(- - l4L-Cfc'i+C+,_Vi+ , /2 = I + +)(- - \4L-^k'L+ct-^q'L+ ■ (20) 

Following the arguments given above, we find that the anomalous exponent for /i is a/2 = 2{2d/7T — 1)^ + 1 > 1, 
while that for /2 is a/2 = 2[28/'k — 1)^ -I- 2 > 1. Therefore, both /i and are irrelevant operators in the charge 
ordered phase. The tunneling does not affect the essential physics (i.e. the quantum phase transition) of the system. 
Its influence is only felt through the phase shift b which depends on the tunneling t. J7kti then, is a function of J7, F, 
and now t. Thus, the tunneling t does affect where the QPT occurs (Fig. 1 in the main text); in particular, to observe 
the QPT in experiments, t can be tuned instead of or in addition to the interaction U' . 
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